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Abstract 

In this paper we analytically study the effects of the lattice discreteness on the electron 
band in the SSH model. We propose a modified version of the TLM model which is derived 
from the SSH model using a continuum approximation. When a soliton is induced in the 
electron-lattice system, the electron scattering states both at the bottom of the valence 
band and the top of the conduction band are attracted to the soliton. This attractive force 
induces weakly localized electronic states at the band edges. Using the modified version 
of the TLM model, we have succeeded in obtaining analytical solutions of the weakly 
localized states and the extended states near the bottom of the valence band and the top 
of the conduction band. This band structure does not modify the order parameters. Our 
result coincides well with numerical simulation works. 
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§1. Introduction 



The trans-polyacetylene is an ideal material which bears a non-linear excitation. A 
soliton is the non-linear excitation which changes the arrangement of double bonds 
of carbons. Su-Schrieffer-Heeger model (SSH model) is the most simple lattice model 
which clarifies the soliton.S' In this paper, we study effects of the lattice discreteness 
on the electron band and the order parameter of the soliton in the SSH model. The 
Takayama-Lin-Liu-Maki model (TLM model) is known as a continuum version of the 
SSH model and gives a analytical solutions for the soliton, the polarono^ and the 
soliton lattice.!^ Takayama et al. derived this model from the SSH model retaining 
only the lowest order term with respect to a/^, where a is the lattice constant in the 
direction of a polymer chain and ^ is the soliton width. To take into account the 
discreteness of the SSH model, it would be a correct way to include higher order terms 
which will be treated in a perturbational way, the TLM model being the unperturbed 
Hamiltonian. Obtained results have to satisfy a selfconsistency condition between the 
order parameters and the single-electron wave functions. 

There were several works which determined an acoustic component of the order 
parameter. For example, Maki obtained the order parameters up to the first order of 
a/^ using the TLM modeLS At low temperature, soliton friction motion is dominated 
by the interaction between the acoustic phonon and the order parameter. Some 
workers used his result to study the soliton diffusion.!^!) Kurita, Ono, and Wada have 
taken into account the fact that the band width is finite, using a model introduced 
by Gammel.&& It retains the property that the electronic wave number is confined 
in the first Brillouin zone. Thus, the finiteness of the electron band reappears in 
a natural way even in the continuum limit. Kurita et al. have obtained the first 
order corrections to the electronic wave functions and the order parameters by the 
perturbation method. They have found that there is a divergence in the second order 
calculations. 

It has been pointed out that the first order acoustic component of the order 
parameter is working as an attractive potential to give rise to electronic localized 
levels at the bottom and the top of the band.&i* There have been several works 
which have solved the SSH model numerically. Kurita et al. have pointed out that 
there are two electronic levels which are localized at the soliton or at the polaron, in 
addition to the well-known mid-gap states.!-* The structure of the electron band for 
a soliton is shown in Fig.l. One of them is below the bottom of the valence band 
and the other above the top of the conduction band. Fu, Shuai, Liu, Sun, and Hicks 
have obtained the two localized levels for the soliton together with two more levels 
in the energy gap.E^P Kurita et al. have analyzed the phase shifts of the electronic 
wave functions of the SSH model numerically in ref. and showed that there are 
two types of electronic states which are classified according to a property roughly 
corresponding to parity. There are accordingly two types of the phase shift. They 
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agree with each other except at band edges. The present author and Wada showed 
in ref . that the weakly locahzed states at the edges of the band can be determined 
analytically by solving equations given by the sum of the TLM Hamiltonian and the 
first order correction terms without using the perturbation method. 

The purpose of this paper is to develop the method to obtain the scattering states 
of electrons near the band edges as well as the localized states analytically up to the 
first order of a/^. We propose a modified version of the continuum model which is 
constructed introducing the electronic wave functions at the edges of the band in 
addition to the electron wave functions at the fermi level. The derived equation for 
the electron wave function contains the acoustic component of the order parameter as 
the attractive potential. The selfconsistency condition is reconstructed using both the 
electron wave functions at the edges of the band and at the fermi level. Solving the 
equation, we obtained analytical expressions of the electronic wave functions. Using 
this model, we will show that the behavior of the phase shift of the extended electronic 
states for a soliton and the existence of the weakly localized states are consistent with 
the small distortion of the acoustic component of the order parameter. 

In §2, the formulation of the modified version of the TLM model is constructed. 
The self-consistency between the lattice order parameter and the electron wave func- 
tion will be shown in §3. We solve the equations in §4 and obtain the structures of 
electronic bands. The phase shifts of the extended states are also derived. 



§2. Continuation of the Su-Schrieffer-Heeger model 

The Hamiltonian of the SSH mo de© IS 

n,s 

+^K^K+i - u^f + ImJ2uI (2.1) 

n n 

where the overlap integral t„+i,n is given by t„,„+i = to — a{un+i — Un)- tn,n+i is a 
linear function of the bond length between the n and n + 1 sites. The quantity ^ is 
the creation operator of a vr electron at site n with spin s, and m„ is the displacement 
of the ra-th site in the polymer chain. The second term is the harmonic potential for 
the bond between the n and n + 1 sites with the spring constant K due to the a 
electrons. The third term is the kinetic energy of the unit, M being the mass of the 
unit. 

The bond variable y{n) is defined by 

y{n) = Un+i - Un- (2.2) 

Considering that the dynamics of lattice is sufficiently slow in comparison to that of 
the electrons at the Fermi level, we have, in the adiabatic limit. 
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Then, the set of electronic wave functions ^k,s{n) is determined as eigen- functions of 
the equations 

Ek^k,s{n) = -{to - ay{n - - 1) - {^o - ay{n)}<^kA^ + 1), (2.3) 

with energy E^. This equation is derived by the variational principle with respect to 
C^. Using ^k,s, the total energy is given by J2'k,s{^k,s\HssH\^k,s) ■ This is minimized 
when y{n) satisfies 

= -;^E {^k/{n)^k,s{n + 1) + c.c} + ^ E E {^k/{n)^kAn + 1) + cc}, 

k,s n k,s 

(2.4) 

where J2 means a summation over the occupied states. The second term comes from 
the constraint 

E?/H = o- (2-5) 



n 



As is shown by Su-Schrieffer-Heeger,B this model bears the bond order wave (BOW) 
state, which has the bond alternation of the wave number n/a. In this state, the 
bond variable y{n) is given by 

Vn = (2.6) 

where A is the optical component of the order parameter. In the SSH model, the 
quantity A is determined by the selfconsistency equation, 

1 = 2A / =d{ka). (2.7) 

•^0 ^(2to)2cos2A;a + A^sin^fca 

When the non-linear excitation is induced in the BOW state, the order parameter 
A is distorted near the center of the excitations. Makii^ had pointed out that the 
acoustic component of the order parameter also suffers a distortion around the non- 
linear excitation. The bond variable is written as follows: 



'M = ^[({na + ^) + {-rK{na + ^)}, (2.8) 



where C, is the acoustic component of the order parameter. Takayama-Lin-Liu-Maki 
have written the electron wave function nearly at the Fermi level by the summation 
over the left-going and the right-going waves, & 

^k,siji'a) = [Ws{na) explikpna} — iVs{na) exp{—ikFna}]y/a, (2.9) 
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with the Fermi wave number kp = Tx/2a. On the other hand, near the bottom or the 
top of the electron band, the electron wave function should be given by 

$fc,,(na) = {Z,M-in(na)e'2'=^("'^)}v/^. (2.10) 
for \E\ ~ 2to, 

where Z and Y are slowly varying functions. Using W , V, Z and Y, the electronic 
states in the presence of the non-linear excitation will be analyzed. Substituting 
eqs.(p.8[) and ( |2.9| ) into eq.(p.3|), and introducing ipk,s{x) by 

^m(^)= f w!? V (2.11) 



Vs{x) 

we obtain the equation for the electron, 

Ek^k,s{^) = -iAoo-39?/'fc,s(a;) + A(x)cri?/^fc,,(x) 

+mil^^{d(ix))asiJk,s + iC(x)a3(9^fc,.(x))| + . . . , (2.12) 

where x = na, and 2Ao being the electron band gap. The quantity mi = a/^ = 
Ao/(2to) is the smallness parameter. The differentiation d is defined hj d = 
where ^ = 2tQa/ is soliton width. In the calculation, ipk,s{x ± a), the functions 
A(x ± a/2) and ({x ± a/2) were expanded around x. 

Substituting eqs.(|2.8|) and ( |2.1(J| ) into eq. (PTB|) and introducing v^5,s(x) by 

^U^)=(i±l), (2.13) 



we obtain the equation for the electron near the edges of electron band, 

Egipg^six) = (^-2to - ^mid"^ + C(x)^(j3V9g,s(a;) 

-imi Q(9A(x)) + A9)aiy.,,,(x) + . . . . (2.14) 

The functions (Pg,s{x ± a), A(a; ± a/2) and ({x ± a/2) were expanded around x. 
Equations( p.ll| ), ( |2.12| ), ( |2.13| ) and ( |2.14| ) are the equations of the electron wave 
function. 

Substituting eqs.( |2.9| ) and ( |2.1CI| ) into eq.(p^), we obtain the selfconsistency rela- 
tion between the electron and the bond variable y{n). Using eq. (^78|) and equating 
the slowly varying term and fast oscillating term by (— )", we obtain the equations 
for the order parameters A (a;) and ({x) 

~ Act^a-r ./ 

k,s 

-mi^-^Y!'^{i^'Pisi^))^i'PiA^) - flA^yii^ipq^six))}, (2.15) 
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+mi^^^'i{(94 3(x))a3V'fc,.(x) - ijl/x)a3{d^k,s{^))}, (2.16) 

where the summation J2' is performed over the occupied states. fq^s{x ± a/2) and 
ipk,s{x ± a/2) were expanded around x. 

Equation (p.l2|) determines the wave function of the electron at the Fermi level, 
while eq.( |2.14| ) is the wave function of the electron near the top of the conduction or 
the bottom of the valence band. Equations ( |2.15|) and ( |2.16|) are the selfconsistency 
conditions for the order parameters A and (. 

The zeroth order terms of the expansion with respect to mi in eqs.( p.l2"[) and 
( |2.15| ) give the TLM model. Equations ( |2.14| ), ( p.l6| ) and the higher order terms 
with respect to mi in eqs. (|2.12|) and ( p.l5|) come from the lattice discreteness of the 
SSH model. 

The order parameter of soliton in the TLM model is given by & 

A(x) = AotanhKx, (2.17) 
with K = 1/^. The electron wave function for the scattering state is 



, , _ ^T, + 7 + itanh(Kx) \ 

^^^^^ 2 U(^7-7-itanh(/€x))y' 

N, = (L(l + f)-l/«:)(-i/2), 



'7 

= ±v/rT7, (2.18) 
7 = efc, (2.19) 

where k is the wave number. In eq.( |2.18| ), the plus sign corresponds to the conduction 
band and the minus sign to the valence band. The energy of the electron is given by 

Ek = AoE^. (2.20) 

The electronic localized state inner side of the band gap is given by 

ipBoix) = ^^sech Kx (^\^ , (2.21) 

Em = 0. (2.22) 

The quantity Aq is determined by the lowest order term of the selfconsistency equa- 
tion ( p.l5| ). This condition becomes 

l = 2Alog(^), (2.23) 
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with A = 4aV(2toi^vr). 

Maki had analyzed the sohton using eqs. (|21^) , (|35D and (|21^),& neglecting the 



terms related to (fq^s and retaining terms up to the first order with respect to mi. He 
obtained the acoustic order parameter 

((x) = -!^Aosech 'ra, +0(mi2). (2.24) 

The optical order parameter and the electronic wave function were unchanged. Nu- 
merical and analytical works in the SSH model have pointed out that there are weakly 
localized states around the soliton in addition to the mid-gap state. Their energy 
levels are located at the top of the conduction band and the bottom of the valence 
band.B&0' The acoustic component of the order parameter C, works as an attractive 
potential at the soliton center, making one weakly localized state at each of the band 
edges.1^ As shown in Fig.l, the level structure is different from that obtained in ref.jQ. 
Therefore, it is important to see whether the selfconsistency between eq.( p.24| ) and 



the electronic states is satisfied or not. In the next section, we obtain the electron 
wave function near the edges of the electron bands and the selfconsistency will be 
confirmed up to the first order with respect to mi. 



§3. Selfconsistency 



In the present section, the electronic states and the lattice configuration are self 
consistently determined, using eqs. (|2.12| ) ~ ( |2.16| ). We expand the electron wave 
function by the smallness parameter mi. 



^kA^) = ^/'i'J(x) + miV'2(a;) + ... 
^g,six) = ^f}ix) + miipl^lix) + ■ ■ ■ 



and the order parameters, in the same way. 



A(x) 
C(x) 



A(°)(x)+miAW(x) + --- 
C^^^^x) + miC^^\x) + ■ ■ ■ . 



Then, the equations for electrons, eqs. ( p. 12]) and ( p.l4[ ), become 









rk,s 








^(0)^(0) 






+ E(0V^'^ 

' q rq,s 





(3.1) 



(3.2) 



(3.3) 
(3.4) 
(3.5) 
(3.6) 
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where the Hamiltonian H^'^\ h^^'' = 0, 1) are given as follows: 

= -Aoia39 + A(o)(x)(Ti, 

rrii 



-(9A(°)(a;)) + A(°)(a;)9 



(Ti. 



(3.7) 
(3.8) 

(3.9) 

(3.10) 



Because the electrons near the edges of the electron bands suffer the attractive po- 
tential given by the acoustic component of the order parameter C{x), the second and 
the fourth term in the right hand side of eg. ( p.9D is included in the unperturbed 
Hamiltonian h^^\ Substituting eqs. (|3.1|) and (|3.2|) into eqs.( 2.15 ) and ( p.l6| ), the 
selfconsistency equations become 



A(°)(x) 
A«(x) 



C«(a;) 



' K 
K 
' K 



«IlVS\a;)aiV^l^°](x), 

k,s 

V)c^iV^S(a;) + V)(Ti^i°](x)) 



(3.11) 



k,s 



q,s 



EVS\a;)a3^W(a;) 



K '^'^ 



Aa' 



(3.12) 
(3.13) 



aE'{^S V)^3^i5(a:) + (a:)a3^(:i(x)} 



q,s 



2a^ 



K 



^E'i{(54f )(^)^3<](x) - V^if (x)a3(94°](x))}. (3.14) 



k,s 



Using eqs.(|3.3|), (|3l5|), ( p.7| ), (p. 91), ( |3.11| ) and ( p.l3| ), an unperturbed solution is 
obtained. 

Equations ( p. 71 ) and ( |3.1lD are the same equations as in the TLM model. Thus, 
in the soliton case, the unperturbed solution ip^^\x) for the scattering state is given 
by eq.( |2.18| ) and the solution for the mid-gap state is given by eq. ( |2.21| ). Our model 
does not modify the electron wave function near the Fermi level which is given by 
the TLM model. 

The optical component of the order parameter is given by 



A(o)| 



X] 



Aq tanh KX. 



(3.15) 
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The band gap Aq is given by the equation 

l = 2Alog(H|^). (3.16) 

where A is introduced as the cut off of the wave number k. 

The wave function of the electron near the bottom of the valence band or the top 
of the conduction band is obtained by eqs. (|3.5|) and (P-Qj). The wave function ip is 
given by 



E, = -2to + eg, (3.18) 

for the valence band, and by 

^r(-) = (J,,). (3.19) 

E, = 2«„-£„ (3.20) 

for the conduction band. The functions (f)q{x) are solutions of the Schrodinger equa- 
tion 

eM^)={^-^Aod' + C^'\x)+m,C^'\x)y,{x), (3.21) 

where C^^\x) and C^^\x) is the real function given by the selfconsistency equation 
(|3.13| ) and (|3.14|) . The quantity Eg is the eigen- value of eq.( p.21| ). As it will be 



shown later, the forms of C^^\x) and C^^\x) do not depend on the structure of 4>q{x). 
Because (pq^x) is the solution of the second order differential equation, there is a 
completeness relation, 

Y.^l{^)Uy) = 5{^-y)- (3.22) 
<? 

where 0g(x) is renormalized as / 0*(x)0g(x)da; = 1. In eq. (|3.22| ), the summation is 
taken over all eigen-functions {(pqi^x)}. From the selfconsistency equation ( |3.13| ) and 
the completeness relation (|3.22|) , the zeroth order part of the acoustic component of 
the order parameter, C*^°^(x), is the constant. Since, from eq. ( p.5| ), the integration 
over X of C,{x) is zero, we obtain 

C^'^\x)=^. (3.23) 

Therefore, substituting eq. ( |3.23| ) into ( p.8| ), the first order Hamiltonian H'^^'> is zero. 
From eq.( 3.4|) , the first order wave function is zero. 

il)^^\x)=Q. (3.24) 



9 



Substituting eqs.( p.l9| ) and ( p.24| ) into eq.( p.l2| ), the first order term with respect to 
mi of the optical component of the order parameter, A^^\x), is given by 

AW(x) = 0. (3.25) 

From eg s . (13.61 ) and ( p.lO ), the first order electron wave function (^^^\x) is given 
by the perturbation. 



—i J dyQq{x,y)(^-sech "^Kx + tanh Kxd^aiip^^\y), (3.26) 



where we used eg. ( p.l5| ). Qg(x,y) is the perturbation kernel, 

= (3^27) 

In the summation X)" over k' , the wave number q = q' where sgn{Eq) = sgn{Eqi) 
is excluded. Substituting egs.( |3.17| ) and ( |3.19| ) into eg.( |3.27| ), the kernel Qq{x,y) for 
the valence band becomes 



Qqix,y) 



1 0\ „ 0g.(x)0;,(y) / \ mi 



5(x-y), (3.28) 



where is the eigenvalue of the Schrodinger eguation ( p.21j ). In the calculation from 
the first line to the second, we approximated + \Eqi\ by 2 x (2to), since eq/{2tQ) 
is the guantity of the order of m\. We used the completeness of (pk{x), eg. (|3.22|) . As 
in the same way as eg.( p.28|) , substitution of eg.( p.l7|) into eg.( p.27|) gives the kernel 
Q{x,y) for the conduction band. 



(3.29) 



,>0 ^Q' 



From egs.( |3.26| ), ( p.28| ) and ( p.29| ), the first order wave function is 



if''q\x) = — i sgn(£',) — f -sech kx + tanh. nxdj ai(f^q\x) . (3.30) 



mi / I 

I t . , . 

2 V2 J 

The perturbation Hamiltonian h^^'^ is expressed by using cti. As we can see from 
egs.(3.17) ~ (3.20), this matrix causes the mixture of the electron wave functions 
at the valence band and the wave functions at the conduction band. Because the 
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energy difference between eqs.(3.18) and (3.20) is the order of 4to, the perturbation 

x)/2. Then, the perturbation 



(0), 



term (^)j> (x) is the order of Aq / (4to) x y?^ 
Hamiltonian h'^^^ gives the second order contribution to the electron wave function 
^q{x). Substituting eqs. (p:2l|) , (|3T7| ) and ( F^ into eq. (|3l4D , C^^\x) is given by 



C 



(1), 



X] 



2 



(3.31) 



while v9^^-'(a;) gives no contribution to the acoustic component of the order parameter 



C{x). This result is the same as eq.( |2.24 ). Because of the completeness of the electron 
wave function ^q{x)^ the electronic states near the edges of the band do not modify 
the order parameter up to the first order with respect to mi = a/^. The order 
parameters are determined by the electrons near the Fermi level. 



§4. Electronic States at the Band Edges 



The electron wave function near the edges of the electron band is given by 
eq. (p.l7|) , where the function 0g(x) is the solution of the Schrodinger equation ( p.21|) . 
Because the acoustic component of the order parameter, in the soliton case, is given 
by eq. (|3.31|) , this equation is written as follows : 



2e 



miAo 



0(x) 



d{Kx) 



sech Kx\(j){x). 



(4.1) 



Using the hypergeometric function, 0(a;) for a positive e is written by 

1 + tanh(Kx) 



05 (x) 

2e 
miAo 



e^^^F[l + /?o,-/?o,l + ieg; 



(4.2) 



where q is the wave number. (3q is defined by eq.( [A.6|) . The details of the calculation 
are given in the Appendix A. From eqs.( p.l7D and ( |4.2| ), the electron energy Eq is 
given by 



±2to( 1 - 



(4.3) 



As eq. ( [4.1|) is symmetric under the transformation x —x, (pq{—x) is also a solution 
with (2£:)/miAo = {^qY- We use combination of 0g(x) and 0q(— x), which are defined 
by 



0g,i(a;) 



1 



L(|l±6(eg)|2 + |a(eg)|2) 



^q{x)±<pq{-x)], 



(4.4) 
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where the upper sign corresponds to i = e (even) and the lower sign to i = o (odd). 
The normalization factor is calculated in the Appendix B. A negative e corresponds 
to the localized solution. It is given by 



r(l + 2/3o) /sech 



A) 



\ er(i + /5o)r(/5o) V 2 



(4.5) 



From eqs.( |3.17|) and 



the energies of the weakly localized states are 

i2 



Ef 



This form is the same as the eq.(4.28) given in ref. | 
the TLM model, the quantity /5o is given by 0* 



/So 



V5 - 7A - 1 



(4.6) 

In the finite band version of 

(4.7) 



In the present theory, the electron cosine band is linearized. Then, the contribution 
of the order A is dropped from (3q and the quantity mi in the present theory is a little 
larger than that obtained in the finite band version of the TLM model. 

Substituting eqs.(|4.4|), (^3|), ( |3.17| ) and ( |3.19| ) into eq.( ^.lU| ), we obtain the elec- 
tron wave function $ 



9' 



^q,i{x) = y/a(f)g^i{x) X 

^/a(j)i{x) X 



9ii ^ for the scattering states 



1, for Eg < 0, 



for the localized states 



(4.^ 



with i = e,o. The wave function $f (x) of the electron at the bottom of the valence 
band is shown in Fig. 2. The band structure is schematically shown in Fig. 3. We have 
selected the cut off of the wave number g as A<^ = kp/2 = 7r/(4a) and the cut off of 
as A = kp/2. 



§5. Phase Shift of the Extended States in the Electron Band 

The localized states corresponding to to those discussed in the previous section were 
also obtained by the numerical calculation of the soliton in the SSH model. It was 
shown that there were two branches in the phase shift of the electron wave function 
in the SSH model.0^ The phase shift for the electron wave function with even parity 
approached — vr as g — 0, while that for the wave function with odd parity approached 
as g — > 0. 



12 



The behavior of the phase shift is reproduced by our analytic wave function, 
eq. ([4.4|) . As \x\ oo, the electron wave function, eq.( [4.4|) , behaves like 



x[aK5)e*l + (6«9)±l)e-"l'l|, (5 A) 

where i = e, o denotes even or odd parity. Dn is defined as 

^ _ I 1, for the valence band , (^2) 

" 1 Q^ikpx^ ^Yie conduction band. 

Equation (|5.1|) is rewritten as 



where 



1 + 

(5.4) 



am - 

1 - bm 



am 

We plot the phase shifts Sg and (5o in Fig.4(a). They are drawn in the region from 
g = to g = A^, where we selected as kp/2. Figure 4(h) shows the phase shifts 
obtained by the numerical calculation in the SSH model.Q) The agreement between 
the analytical result and the numerical one is remarkablly good. The 6e approaches 
— TT as g — >■ 0, while So approaches as g ^ 0. The both phase shifts, 6e and 6o increase 
and approach to at the cut off g ~ A^. These behaviors indicate that of the electron 
plane waves near the edges of the electron band are pulled into the soliton center by 
the attractive potential ({x). This result is consistent with Levinson's theorem. The 
theorem says that the number of localized states at the bottom of the band reflects 
the phase shift. If the numbers of localized states with even and odd parities are 
rif, and respectively, the phase shifts of the even and odd extended states go to 
— (2ne — l)vr and — 2no7r, respectively, as g ^ 0. 0^ The present case corresponds to 
TT-e = 1 and Ho = 0. 

The phase shift near the Fermi level is given by 

e'' = (5.5) 
1 + 7 
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where 7 is the quantum number defined by eq.( |2.19| ). This phase shift ( p.5| ) is the 
same as that obtained in the TLM modeL We have to be careful about the definition 
of the wave number. At the Fermi level, k{= 7/O is zero. Thus k is considered as 
the deviation of the wave number from kp. Therefore, 

k = q — kp- (5.6) 

The phase shift 5 in the valence band is plotted as the function of q in Fig. 4(a). 

= corresponds to g = kp and the cut off A; = —A = —kpjl corresponds to 
q = kp — A = kp/2, respectively. The phase shift 6 increases as q becomes larger, 
approaching to 71/2 as q — >■ kp. This behavior is also in good agreement with the 
numerical result shown in Fig.4(b). There is a discontinuity of the phase shift at 
q = A^. This is because we have used eq.(2.11) for < g < A,^ and eq.(2.13) for 
<q < kp. 

§6. Summary and Discussion 

We have studied the effects of lattice discreteness on a soliton in the SSH model 
using a modified version of the continuum (TLM) model. We solved this model in 
the adiabatic limit. 

When there is a soliton in the SSH model, the arrangement of the double bonds 
is changed and the ions near the soliton are shifted to the soliton center. Then, the 
acoustic component of the order parameter, C{x), has non-zero amplitude near the 
soliton center and it acts as a potential for the electrons at the band edges. Since 
({x) acts as an attractive potential, weakly localized states appear at the band edges. 
The levels are very close to the corresponding band edges. Because the amplitude 
of the acoustic component is very small (of the order koAq), the states are more 
weakly localized than the mid-gap state. In the TLM model, this effect (the acoustic 
component) is not included because the interaction between the acoustic component 
({x) and the electrons at the bottom of the valence band is omitted in the continuum 
limit. 

To improve the TLM model, we have taken into account the electronic wave func- 
tion at the bottom of the valence band and the top of the conduction band in addition 
to the electron wave function at the Fermi level. We have considered the discreteness 
by taking into account the terms of higher order in a/^, when the expansion parame- 
ter a/C, is sufficiently small. We have studied the electronic states of electrons, using 
the unperturbed Hamiltonian eqs. (|3.7|) and (|3.9| ), which are composed of the zeroth 
and a part of the first order terms with respect to mi. The selfconsistency between 
the electrons and the order parameters have been satisfied up to the first order of 
mi. 

The wave functions of weakly localized states are given by eq. ( [4.8|) . The energy 
levels of the localized states are given by eq. ( [4.6| ). This result has the same form as 
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given in the finite band version of the TLM modehl^ Concerning the scattering states, 



we have obtained the dispersion relation. eqs. (|2.2CI|) and (^4.3|) , which corresponds to 
the scattering states in the SSH mo 



Eq = ±V (2to)2 cos2(ga) + sin'^iqa). (6.1) 



The wave functions are given by eqs.( [4.2| ) and ( [4.4| ). Using these wave functions, the 
phase shift analysis has been performed for the two groups of the electronic states, 
with even and odd parities respectively. The two phase shifts coincide with each other 
over a wide range of the wave number. In the regions of small wave numbers, the 
two take different values which are consistent with the fact that there are localized 
states with even parity. 

Using the present model, the order parameters, eqs.( p.l7D and ( p.24|) , which were 



given by ref . |^ , have been reproduced. In the expression of the electron wave func- 
tions near the band edges, eq. (|4.^ ), the complete set of the function {(pq,i{x)} have 
been introduced . Because of the completeness of {(pg^i^x)}, the acoustic component 
of the order parameter, C{x), is unmodified by the weakly localized states at the edges 
of the band. The order parameters are determined by the electrons near the Fermi 
level. 

We have made a modified version of the TLM model and solved it perturbationally 
up to the first order with respect to a/^. The present model has a merit that it can 
deal with the electrons near the bottom of the band as well as those near the Fermi 
level. It has another advantage that we can treat the interaction between the electron 
and the acoustic component of the order parameter. Thus, this model is suitable for 
analyzing the interaction between the acoustic phonon and the non-linear excitation. 
It is interesting to study the diffusion motion of the nonlinear excitation induced by 
the interaction using this model. 

To obtain the exact result up to the second order of K,a, we should make a higher 
order perturbation calculation. However, since the cosine band in the SSH model is 
linearized by the continuation procedure, the contributions of the order A ■ mi, ml 
and higher which come from the scattering states are omitted. Kurita et al. obtained 
the acoustic component, C{^)i using the finite band version of the TLM model.!) C(x) 
was given by 

C{x) mi - ^A^ Ao cosh~^ kx + O(m^). (6.2) 

Then, our result of the acoustic component, eq.( |3.31|) , is a bit larger than that of the 
SSH model. To make higher order perturbation selfconsistently, we should construct 
the unperturbed Hamiltonian using the finite band version of the TLM model.i''§-'0 
It is a bit tedious but important problem to explore. 
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APPENDIX 



Appendix A: Solutions of the Differential Equation (4.2) 

In this Appendix, we will consider the differential equation 

d^ 

e'(f)(x) = —— — r:70(a;) — sech (A.l) 

First, the localized solution with a negative e' will be considered. Introducing t{6) 
and 6 by 

= (1 - tanh^ Kx)^t{e), (A.2) 

e = 1±I|^. (A.3) 

We write eq. (|A.l|) in the form 



^(1 - 0)^ + (1 + v^)(l - 2^)^ - i2R - 5' + V^)\m = 0. (A4) 

This is the hypergeometric differential equation.El'' With the help of eqs.( |A.2| ) and 
(|A.3| ), the solution ^^(x) is given by a hypergeometric function 



M^) = mi-e)]^F[V^-/3o,V^ + (3o + l,V^+l,0], (A.5) 

where 

(3o = {V5-l)/2. (A.6) 

As we can see from eq. ( |A.5| ), (f)e{x) approaches zero, as x goes to —oo. In order that 
4>e{x) is finite in the limit of x — > oo, the hypergeometric function F should be a 
polynomial. So we can put \/—e' — (3o = —n with a non-negative integer n. Since /?o 
is about 0.6, n should be zero. Then we have 



V^ = Po. (A.7) 
As Po is positive, we have a localized solution 

Mx) = [eil-e)]^F[0,2(3o + l,(3o + l,0] 



/sech Kx\ 



/3o 



) ■ (A.8) 



V 2 

The eigen value e' is given by 

e' = -0112. (A.9) 
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To obtain the solutions of scattering states with a positive e', we replace \/—e' in 
eq. (|A.5|) by 17, that is e' = 7^, to get 

= [^^(1 - e)pF[-i3^ + i^, + 1 + i^, 1 + e]. (A.io) 

The hypergeometric function has the following relations 

F[a,/3,7;^] = (1 - /^F^ - a, 7 - A 7; (A.ll) 



+ /3, a + /5 - 7 + 1; 1 - 



r(7-a)r(7-/5)' 



7— /3f[^ _ c,,^ _ _ a _ + 1; 1 _ ^] 

(A.12) 



We write (f)^ as a function of 1 — 6*, using eqs.( A.lI ) and ( [A.12 ) and the relation 



= to get 



= e^-^^F[l + /3o,-/3o,l + i7;^] 



(A.13) 



where 

0(7) 
6(7) 

We rewrite ^^(a;) by g = 7/^, 



r(i + i7)r(i7) 



r(l+/?o + i7)r(-/3o + i7)' 

r(i + i7)r(-i7) 



(A.14) 



(A.15) 



q is the wave number of the electron at the band edges. From eqs.( [A.13D and ( [A.15| ), 
we can see that the function (p^ has the asymptotic forms 



as X ^ —00, 



a{^q)e^''^ + b{^q)e '^"^ as x 00. 



(A.16) 



Appendix B: Calculation of Normalization Factor 

Consider an hermitian operator Lq and its eigen function u^{x) 



Lo 



(B.i; 
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where the eigenvalue -£^(7) is related to the quantum number 7 by £'(7) = 7^. We 
differntiate eq. (p.l| ) by 7 and multiply by u*_^ from the left to obtain 

M;(Lo-7'K = 27w>r (B-2) 
The hermitian conjugate expression of eq. (p3.1|) gives 



{Lqu* - 'y^u*)u = 0. 



The difference between the two is integrated to give 

= 27 /" utfU^dz. 

J b 



'd /\ 



Suppose we put 

which give (f)q^i{x). With the help of eq.( |A.16| ), we get 



(B.3) 

(B.4) 
(B.5) 



'x) Pdx 



(|l±6(7)r + |a(7)r)L + 0(l) 



(B.6) 



The wave functions of the localized states eq.( |A.5| ) can be obtained from u^{x), if we 
perform an analytic continuation 



7 ^ 7i = 2(a - n), 
where is a non-negative integer. Equation (|B.4|) becomes 



271 / {u^iYdz. 



'dz""^ ^'^ dz^^^' 

Since eq.( |A.14| ) gives 

a(7i) = 0, 

6(71) = i(-l)", for 71 = i(3o 
substitution of eq.( |A.16|) with 7 = 71 into eq. (|B.8|) gives 

iV-2 = ^'^^\cP^,^^{x)Ydx 

J-L/2 

= ia'(7i)^(7i)/« 

n\r{(3o + l-n)r{(3o-n) 



:b.7) 



(B.8) 



r(l + 2(3o - n) 



In particular, for 



L/2 
-L/2 



COS 



h(nx)-^^^dx - 2^'°^' r(/?o + i)r(/?o) 



(B.9) 



(B.IO) 
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Figure Captions 



Fig. 1: Levels of localized states of a soliton. 

Two states, B and C, are close to the band edges and weakly localized. 
The energy of the these two states are given by eg. (^^) . The state A is 
a mid-gap state. 

Fig. 2: Electron density distributions of the bound state wave functions for the 
coupling constant A = 0.2. 

The ordinate is scaled by 1/.^. The abscissa is x/^. is the electron 

density at the bottom of the valence band. |$mP is the electron density 
in the mid-gap localized state 

Fig. 3: Schematic electronic band structure of a soliton. 

The dispersion relation at the Fermi level, eq. (|2.2(]| ), is plotted as the 



function of /c = q — kp, which is identical with that of the TLM model. 
The dipsersion relation at the band edges, eq.( |4.3|) , have a parabolic form, 
which is plotted as the function of q. We selected the cut off as A = = 
kF/2. 

Fig. 4: Phase shifts of the wave functions as a functions of q. The dimensionless 
coupling constant A is selected as 0.2. The ordinate is scaled by l/vr. The 
abscissa is q/kp- 

a) The phase shifts 6e and 6o given by eq.( |5.4|) . 6 is the phase shift of 
electrons near the fermi level. We selected the cut off as A = A,^ = kF/2. 

b) The phase shift of electrons in the SSH model obtained in ref. 0. 
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